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THE  ZENITH  ANGLE  DEPENDENCE  OF  COSMIC  RAY  PROTONS 


CHAPTER  I 
INTRODUCTION 

The  composition  and  spectral  distribution  of  the 
cosmic  radiation  is  a function  of  atmospheric  depth,  mag- 
netic latitude  and  the  direction  of  observation.  Here  one 
refers,  of  course,  to  the  total  radiation.  A complete 
description  of  the  radiation  would  give  the  intensity  and 
spectral  composition  of  each  class  of  particle  as  a function 
of  these  parameters  and  the  time.  Such  a complete  descrip- 
tion is  not  possible  with  present  experimental  techniques. 

In  successively  refined  experiments  (1-7)  conducted  by  this 
laboratory,  a technique  of  proton-meson  analysis  of  the  radi- 
ation has  been  developed  and  the  spectral  composition  of 
these  radiations  studied  at  sea  level  and  at  3.4  kilometers 
altitude.  These  studies  were  all  for  radiation  of  vertical 
Incidence  at  the  point  of  observation.  It  then  became 
desirable  to  undertake  further  observations  in  order  to 
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find  the  change  in  composition  and  intensity  with  direction 
of  observation  at  fixed  atmospheric  depth.  Two  programs  of 
this  kind  have  been  carried  out.  In  the  work  of  Charbonnier 
(6),  the  spectral  composition  of  the  radiation  of  range 
greater  than  15  cm.  of  lead  was  studied  as  a function  of 
zenith  angle.  As  was  shown,  the  radiation  defined  in  this 
way  consists  almost  entirely  of  mu  mesons.  As  a result  of 
these  measurements,  it  was  possible  to  give  a much  more 
detailed  experimental  check  of  the  theory  of  atmospheric 
development  of  the  mu  meson  radiation  than  had  been  pre- 
viously available.  The  present  work  reports  a similar  study 
of  the  radiation  of  range  less  than  15  cm.  of  lead.  This 
radiation  contains,  in  addition  to  low  energy  mesons,  a 
large  number  of  protons.  Precautions  taken  to  exclude  the 
soft  component  are  described  in  Chapter  II,  Section  1.  In 
the  past,  of  course,  many  studies  of  the  zenith  angle  depen- 
dence of  the  total  hard  and  soft  components  have  been  made. 
The  present  work  Is  the  first  study  of  its  kind  for  any  part 
of  the  nucleonic  radiation.  As  protons  and  mesons  at  lower 
altitudes  are  not  thought  to  be  closely  related  genetically, 
it  was  expected  that  quite  different  dependences  on  zenith 
angle  would  be  found  for  the  two  types  of  radiation.  As  it 
turned  out,  the  angular  dependences  of  protons  and  mesons 
were  found  to  be  the  same,  within  experimental  error.  As 
will  be  shown,  this  must  be  considered  a fortuitous  result 
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of  the  particular  conditions  of-  observation. 

Charbonnier‘8  observations  on  mesons  were  easily 
interpreted  because  of  the  almost  complete  lack  of  nuclear 
interaction  of  these  particles.  Their  distribution  as  a 
function  of  depth  and  angle  can  be  interpreted  entirely  in 
terms  of  the  well  understood  processes  of  decay  and  energy 
loss  by  ionization.  Such  is  not  the  case  for  protons  or 
other  particles  showing  strong  nuclear  interaction.  Cross- 
sections  are  such  that  on  the  average  a dozen  or  more  colli- 
sions will  have  taken  place  between  the  top  of  the  atmosphere 
and  the  point  of  observation  of  3.4  kilometers.  These  colli- 
sions involve  pi  meson  production  along  with  production  of 
secondary  nucleons  at  large  angles  to  the  direction  of  the 
primary.  The  development  of  this  radiation  involves  a cas- 
cade process  qualitatively  similar  to  that  of  the  soft  com- 
ponent. This  has  been  discussed  by  Messel,  Janossy  and 
others  (8-20) . However,  almost  complete  lack  of  either 
experimental  or  theoretical  knowledge  of  cross-sections  for 
the  many  processes  involved  causes  the  result  to  be  one  that 
can  hardly  be  expected  to  be  experimentally  confirmed.  In 
addition,  a simplified  one  dimensional  cascade  has  been 
considered  so  that  the  result  cannot  be  compared  with  the 
present  studies  of  angular  dependence.  It  had  been  expected, 
because  of  the  Isotropy  of  the  higher  energy  proton  radiation 
plus  the  large  angle  scattering  in  successive  nuclear  collisions. 
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that  i he  observed  proton  radiation  at  lower  altitudes  also 
would  oe  almost  isotropic.,  It  was  surprising  to  find  that 
the  intensity  decreased  by  a factor  of  3-2  In  going  from  the 
vertical  to  a zenith  angle  of  45  degrees.  This  result  in 
effect  places  an  upper  limit  on  the  average  deviation  suf- 
fered in  a nuclear  collision.  Rather  than  attempt  a 
detailed  cascade  calculation  based  on  questionable  assump- 
tions regarding  cross -sections , a simplified  phenomenological 
calculation  has  been  made  to  give  a description  of  the  observed 
proton  zenith  angle  dependence. 

As  a by-product  of  these  observations,  additional 
information  on  the  proton  spectrum  has  been  obtained,  as 
well  as  data  on  the  meson  zenith  angle  dependence  for  energies 
lower  than  those  observed  by  Charbonnier. 


CHAPTER  II 


THE  EXPERIMENT 

1 . Experimental  Arrangement 

A sketch  of  the  experimental  arrangement  is  shown  in 
Figure  I . The  entire  assembly  may  be  rotated  about  a North- 
South  horizontal  axis,  perpendicular  to  the  plane  of  the 
figure,  so  that  the  cosmic  ray  beam  at  any  zenith  angle  up 
to  50°  may  be  accepted.  Anti-coincidence  counters  Al,  which 
lie  outside  the  acceptance  beam  defined  by  counters  Cl  and 
C2,  serve  to  limit  the  number  of  cases  in  which  shower  events 
are  recorded.  Only  events  of  the  type  Cl  + C2  - Al  - A2  - 
A3  are  recorded.  A one -centimeter  thick  absorber  In  the  form 
of  a half  cylinder  of  lead  placed  above  counter  C2,  plus  the 
absorber  represented  by  the  walls  of  the  cloud  chamber  and 
the  coincidence  counters  establishes  the  minimum  range  and 
momentum  of  accepted  particles.  The  15  centimeters  of  absorber 
between  C2  and  A2  and  3»  of  course,  determines  the  maximum 
range  accepted.  As  discussed  by  Todd  (5),  such  range  defin- 
ition serves  to  resolve  protons  and  mesons  into  distinct 
momentum  bands,  in  the  absence  of  scattering  and  nuclear  inter- 
action. To  the  extent  that  these  effects  are  ignored,  then, 
the  apparatus  accepts  mesons  of  momentum  0.09  to  0.34  Bev/c, 
and  orotons  in  the  momentum  band  0.4  to  1.0  Bev/c.  In  other 
respects  the  apparatus  is  the  same  as  that  used  by  Garrison 
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(7) . The  inclusion  of  the  one  cm.  absorber  over  C2  has 
removed  a small  overlap  between  the  maximum  meson  momentum 
accepted  and  the  minimum  momentum  of  proton  acceptance 
present  in  his  experiment.  This  modification  is  not  so 
large  as  to  prevent  comparison  of  present  with  previous 
work. 

As  shown  in  Figure  I,  anti-coincidence  coverage 
below  the  15  cm.  of  lead  is  quite  complete.  A double  layer 
of  counters  has  been  used  below  the  absorber  to  prevent 
particle  leakage  between  the  counters.  Along  the  sides,  a 
single  layer  was  considered  adequate,  as  in  any  case  few 
particles  will  suffer  such  wide-angle  scatterings.  To 
prevent  amplifier  overloading  from  too  high  a counting 
rate,  three  anti-coincidence  input  amplifiers  have  been 
provided.  One  was  fed  from  counters  A1  while  the  tubes 
below  the  15  cm.  absorber  were  divided  into  two  groups 
feeding  the  other  two  channels.  The  counter  groups  so 
defined  could  be  separately  scaled  to  detect  defective 
counters . 

As  in  previous  work,  the  2.5  cm.  of  lead  absorber 
above  A1  serve  to  insure  recognition  of  electron  showers. 
Such  showers  either  actuate  one  of  the  counters  Al,  in 
which  case  the  event  is  not  recorded,  or  are  recognized  by 
their  appearance  in  the  chamber.  This  will  be  discussed 


later. 
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2.  Experiments  Performed 

Three  sets  of  measurements  were  taken  with  this 
apparatus  located  at  Climax*  Colorado,  elevation  3*3  km.  and 
magnetic  latitude  48°  North  and  in  continuous  24-hour  opera- 
tion during  the  months  of  August,  September,  and  October, 

1951*  Initial  conjectures  as  to  the  results  of  the  zenith 
angle  measurements  led  to  the  conclusion  that  appreciable 
differences  in  Intensity  would  require  rather  large  angles 
of  observation,  so  for  this  reason,  an  angle  of  45  was 
chosen.  The  three  sets  differed  only  In  the  zenith  angle 
to  which  the  apparatus  was  tilted,  a measurement  being  made 
at  vertical  incidence,  at  45°  to  the  East  and  45°  to  the 
West  of  the  zenith.  Since  the  magnet  could  be  tilted  in 
only  one  direction  with  respect  to  the  trailer  in  which  It 
was  installed,  the  trailer  was  oriented  to  give  an  East  or 
West  zenith  angle.  In  each  case,  a reference  structure  on 
the  magnet  was  made  level  in  order  to  obtain  the  same  amount 
of  tilt  East  or  West.  The  direction  of  tilt  was  changed 
periodically  so  that  all  data  in  a given  direction  were  not 
obtained  at  once . Time  did  not  permit  the  taking  of  data 
at  more  than  one  angle  to  the  East  and  West  plus  the  vertical. 

For  the  vertical  direction,  a total  of  3482  pictures 
was  taken,  of  which  1387  were  suitable  for  measurement  accord- 
ing to  the  selection  criteria  to  be  described  later.  The 
West  measurements  provided  1126  measurable  tracks  out  of  a 


total  of  3962  pictures,  the  East  1304  out  of  4488,  for  a 
total  of  11,932  pictures.  The  magnetic  field  was  reversed 
occasionally  to  indicate  any  asymmetry  from  this  source.  The 
data  for  vertical  and  East  zenith  angles  are  divided  about 
equally  between  field  direct  and  field  reverse,  but  unfor- 
tunately the  West  data  were  all  obtained  with  the  field 
direct. 

3*  The  Equipment 

A brief  description  of  the  principal  components  of 
the  equipment  is  Included  for  the  sake  of  completeness  only, 
as  their  detailed  description  has  appeared  before  (4, 5, 6, 7). 
It  should  be  mentioned  that  the  use  and  constant  improvement 
of  this  equipment  by  the  above  investigators  during  the  three 
years  preceding  the  date  of  this  experiment  was  responsible 
for  smooth,  continuous,  and  dependable  operation  throughout 
the  experiment . 

Figure  I shows  the  essential  features  of  the  magnet 
assembly.  The  cloud  chamber  proper  is  centered  between  the 
magnet  poles,  a conical  hole  in  the  forward  pole  providing 
access  for  photography,  and  a hole  in  the  rear  pole  providing 
for  a tube  connecting  the  sensitive  volume  of  the  chamber 
with  the  expansion  diaphragm  at  the  rear.  Illumination  of 
the  ohamber  was  provided  by  four  Sylvania  type  R4330  flash 
lamps,  whose  light  was  focused  by  cylindrical  lenses  into 
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the  chamber,  and  by  two  300-watt  Incandescent  lamps  for 
visual  Inspection.  These  accessories,  plus  the  necessary 
Geiger  counters  and  absorbers,  were  bolted  to  the  inner  face 
of  the  rear  magnet . The  gap  of  2i  inches  between  the  inner 
magnet  faces,  where  the  sensitive  volume  of  the  chamber, 
Geiger  tubes,  etc.  were  located,  was  maintained  by  steel 
spacers  at  the  right  and  left  sides  of  the  magnet.  These 
spacers  also  served  to  complete  the  magnetic  circuit. 

The  chamber  was  filled  with  enough  60-40  n-propyl 
alcohol-water  mixture  to  saturate  the  volume,  and  finally 
filled  with  argon  to  a total  pressure  of  20  pounds  per 
square  inch  above  atmospheric.  Expansion  is  accomplished 
by  releasing  compressed  air  from  behind  a neoprene  diaphragm 
through  two  pop  valves  electrically  connected  in  parallel. 

The  two  pop  valves  provide  faster  air  release  than  the  single 
pop  valve  used  previously,  while  the  electrically  parallel 
connection  insures  that  the  two  valves  open  simultaneously. 

The  magnets  provide  a field  of  8200  gauss  with  an 
energizing  current  of  800  amperes,  and  are  cooled  by  oil 
flowing  through  the  windings  at  a rate  of  60  gallons  per 
minute.  The  oil  in  turn  is  cooled  in  an  oil-water  heat 
exchanger.  The  cooling  water  was  obtained  from  an  artificial 
lake  whose  constant  head  and  temperature  provided  ideal 
chamber  operating  conditions.  Further  temperature  stabili- 
zation was  provided  by  a valve  controlling  a bypass  on  the 
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heat  exchanger.  Temperature  drifts  were  detected  by  the 
unbalance  In  an  AC  bridge,  one  leg  of  which  was  immersed  In 
the  oil  system.  The  error  signal  thus  derived  was  then 
applied,  through  suitable  circuits,  to  the  task  of  opening 
and  closing  this  valve.  This  system  assured  temperature 
stability  to  0.01°  Centigrade. 

Appropriate  electronic  equipment  provided  for  the 
following  functions. 

a.  Event  Selection.  Simple  twofold  coincidence  of 
resolving  time  4 microseconds  is  obtained  using  a Rossi  type 
coincidence  circuit.  Its  output  is  delayed  by  a 5 micro- 
second univibrator  whose  output  feeds  a second  Rossi  circuit 
in  which  anti-coincidence  mixing  occurs.  An  anti-coincidence 
pulse  of  ten  microseconds  duration  then  arrives  at  the  mixer 
five  microseconds  ahead  of  the  coincidence  pulse,  thus 
blocking  its  further  progress.  The  anti-coincidence  circuit 
has  the  desirable  feature  of  being  continuously  retriggerable . 
Thus,  anti-coincidence  blanking  persists  for  ten  microseconds 
after  the  last  pulse  received,  no  matter  how  frequently  they 
arrive . 

b.  Programming.  This  unit  contains  power  supplies 
including  Geiger  voltage,  a scale  of  16,  with  which  each 
channel  of  the  event  selector,  as  well  as  the  output  of  the 
anti-coincidence  train,  can  be  separately  scaled,  and  mechan- 
ical registers  for  recording  the  C and  C-A  events. 
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c.  Cloud  Chamber  Control.  The  desired  event,  selected 
by  the  event  selector,  and  recorded  by  the  program  unit,  then 
triggers  this  unit,  which  carries  out  the  following  sequence 
of  operations: 

1.  Expands  the  chamber 

2.  Removes  the  clearing  field  voltage 

3.  Records  the  expansion 

4.  After  a time  delay  of  0.06  seconds  to  allow  time 

for  droplet  growth,  fires  the  flash  lamps 

5.  Recompresses  the  chamber  and  reapplies  the  clearing 
field  voltage 

6.  Provides  circuit  sterilization  for  a period  of  90 
seconds  to  allow  time  for  the  chamber  to  recover 

7.  Move 8 the  film  to  the  next  frame 

8.  Monitors  all  supply  voltages 

9.  Provides  for  manual  operation  of  the  chamber  and 
camera  rewind,  and  sterilization  to  avoid  the 
recording  of  spurious  data,  e.g.,  in  the  event 
of  failure  of  the  magnetic  field. 


CHAPTER  III 


CURVATURE  MEASUREMENTS  AND  THE  DETERMINATION  OP 

ABSOLUTE  RATES 


The  determination  of  momentum  of  a particle  is  accom- 
plished by  measuring  the  radius  of  curvature  of  the  track 
and  applying  the  formula ; 


p - 3Br  x 10" 


where  p is  momentum  in  Bev/c,  B is  the  magnetic  induction  in 
gauss,  and  r is  the  radius  of  curvature  in  meters.  With  a 
magnetic  field  of  8200  gauss  as  used,  p « r/4.  Curvature 
measurements  of  nearly  4000  tracks  required  a compromise 
between  inaccuracy  and  the  consumption  of  time.  The  method 
used  involves  the  comparison  of  the  curvature  of  a track 
whose  image  is  projected  on  the  screen  of  a comparator,  with 
that  of  a standard  arc . The  process  whereby  these  arcs  were 
prepared  and  a complete  description  of  the  measurement  tech- 
nique have  been  presented  in  previous  work  (4,5).  More 
recently,  the  optical  system  of  the  comparator  underwent  an 
overhaul,  and  the  magnetic  field  reversals  mentioned  in 
Chapter  II,  section  2,  were  performed  to  disclose  any  curva- 
ture distortion  Introduced  by  the  overhaul.  The  data  pre- 
sented here,  consisting  only  of  low  momentum  particles. 
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showed  no  such  effect.  However,  the  succeeding  experiment 
by  Charbonnier  (6),  In  which  high  energy  mesons  were  accepted, 
was  far  more  sensitive  to  such  an  error  and  did  indeed  reveal 
the  presence  of  one.  Both  the  source  and  correction  of  this 
fault  have  been  explained  exhaustively  by  Charbonnier.  The 
error  amounts  to  introduction  of  an  additional  spurious  curv- 
ature of  1/40  reciprocal  meters  in  such  a way  as  to  increase 
the  curvature  of  tracks  appearing  convex  upwards  in  the  com- 
parator. As  described  by  the  above  authors,  the  preparation 
of  the  standard  arcs  precisely  cancels  out  the  error  for 
tracks  of  this  particular  kind.  For  tracks  of  the  opposite 
sign  of  curvature,  however,  the  compensating  error  inherent 
in  the  standard  arcs  no  longer  compensates,  but  aatually 
accentuates  the  discrepancy  to  the  extent  of  l/20  reciprocal 
meters,  i.e.,  twice  the  distortion.  Hence,  the  tracks  appear- 
ing convex  downwards  in  the  comparator,  and  th®se  include 
negative  particles  with- the  magnetic  field  direct  and  positive 
particles  with  the  field  reversed,  must  have  their  measure- 
ments corrected  by  this  amount.  In  Charbonnier ' s case,  this 
correction  was  necessary.  In  this  case,  it  is  relatively 
unimportant,  but  it  is  felt  that  the  correction  should  be 
made  because  of  the  systematic  nature  of  the  error. 

The  apparatus  used  in  this  experiment  does  not  record 
absolute  rates.  They  can  be  Inferred,  however,  provided  we 
make  several  assumptions.  An  accurate  record  was  kept  of  the 
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so-called  C-A  rate,  which  is  the  rate  at  which  counters  Cl 
and  C2  provide  simultaneous  outputs  (within  4 microseconds) 
unaccompanied  by  anti-coincidence  counts.  The  cloud  chamber 
does  not  respond  to  C A counts  which  occur  during  the  90 
second  sterilization  period  provided  to  allow  thermal  recovery 
of  the  chamber  after  each  expansion.  Neither  does  each  C-A 
count  represent  a single  particle  traversing  the  telescope, 
since  shower  protection  is  not  complete.  Identity  of  the 
nature  of  the  event  causing  the  count  is  learned  on  examining 
the  photographs,  which  are  placed  in  one  of  the  following 
three  categories. 

a.  Single  Tracks,  Such  a photograph  is  characterized 
by  the  appearance  of  a single  track  in  such  a direction  as  to 
indicate  traversal  of  both  C counters.  The  length  of  the 
track,  as  well  as  the  occurrence  of  obviously  pre-  or  post- 
expansion  tracks  Is  ignored. 

b.  Showers.  These  are  photographs  showing  two  or  more 
time-coincident  tracks. 

c.  Misses.  Occasionally  a photograph  will  show  a com- 
pletely empty  chamber.  This  Is  also  interpreted  as  a shower 
phenomenon . 

It  is  assumed  that  the  C-A  counts  occurring  during 
the  90-second  sterilization  time  represent  the  same  statis- 
tical distribution  of  events  as  those  photographed.  This 
means  that  the  single  track  rate  can  be  computed  from: 
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Single  Track  Rate  - (c-A  Rat.) 

This  single  track  rate  then  represents  the  total  rate  cor- 
rected for  the  90-second  recovery  time  and  the  occurrence  of 
showers.  The  statistical  uncertainty  in  this  rate,  as  quoted 
in  the  data  sheets,  is  the  sum  of  the  uncertainty  in  the  C-A 
rate  and  the  statistical  uncertainty  in  the  number  of  single 
tracks.  However,  not  all  single  tracks  are  measured,  since 
we  have  assumed  that  tracks  whose  length  does  not  exceed  70# 
of  the  chamber  diameter  cannot  be  measured  accurately.  As 
compared  with  straight  tracks,  this  criterion  discriminates 
against  curved  tracks  to  the  extent  of  about  3#  at  the  lowest 
momentum  measured  (0.09  Bev/c) . Above  .25  Bev/c  the  discrim- 
ination is  negligible.  Thus  we  compute. 


Number  of  Single  Tracks  of  Momentum  p in  dp  per  unit  time 

(Number  of  Measured  Tracks  of  mom.  p in  dp)  fOA ^ 

- 1 (Total-Hunger  of  Measured  Tracfe) (Single  Track  Rate) 

If  we  divide  this  quantity  by  the  product  of  the  area  of  and 
solid  angle  subtended  by  the  telescope,  herein  called  the 
aperture,  and  the  width  of  the  momentum  interval,  dp,  we 

p 

arrive  at  the  absolute  rate  in  particles  per  Bev/c  cm  sterad 
hour.  If  we  use  the  formula  on  page  13,  p - r/4,  then  dp  - 
dr/4.  For  convenience,  a normalization  factor  can  be  defined 
for  each  set  of  data: 

Normalization  Factor  - {Total  , 

(Aperture) 
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Thus: 

(No.  of  Tracks  of  Mom.  p In  dp) 
Absolute  intensity  - * 

The  only  ill-defined  quantity  included  above  is  the 

aperture  of  the  telescope.  Its  uncertainty  traces  to  the 

uncertainty  in  the  sensitive  volume  of  the  Qeiger  tubes. 

From  the  size  and  separation  of  the  counters,  this  coverage 

2 

was  estimated  to  be  about  1 cm  steradian,  a figure  later 
corrected  to  0.94  (6)  by  comparison  with  the  very  accurate 
determination  of  the  vertical  intensity  made  by  Qreisen  as 
quoted  by  Rossi  (25) . 


CHAPTER  IV 


DATA 

Analysis  of  the  data  from  the  three  sets  of  measure- 
ments appears  in  Tables  I - III.  The  direction  of  the  mag- 
netic field  was  reversed  occasionally  for  the  purpose  of 
detecting  the  presence  of  any  systematic  curvature  distortion. 
As  mentioned  in  Chapter  III,  such  a distortion  is  present, 
necessitating  separate  treatment  for  the  two  kinds  of  data, 
field  direct  and  field  reverse.  In  the  45°  West  case,  only 
one  direction  of  field  was  used,  so  that  here  no  separate 
treatment  is  necessary. 

Listed  In  the  tables  are  the  total  number  of  frames 
inspected  and  the  number  of  frames  containing  showers,  misses 
and  single  tracks  respectively.  -The  number  of  frames  con- 
taining single  tracks  whose  length  exceeds  the  lower  limit  of 
10%  of  the  chamber  diameter  is  listed  as  ’’Total  Measured." 

The  latter  are  further  broken  down  into  positive,  negative, 
and  "infinite,  no  sign,"  these  latter  being  in  the  measurable 
class,  i.e.,  long  enough,  but  of  such  small  curvature  that  no 
definite  assignment  of  sign  could  be  made. 

In  one  case  (45°  East,  field  direct),  chamber  operation 
was  so  poor  for  a short  period  that  accurate  measurements  were 
impossible.  However,  the  tracks  photographed  during  this 
short  time  were  Inspected  and  classified  as  to  type  of  event 
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(single  track,  shower  or  miss,  positive  or  negative).  Thus, 
the  figures  for  these  classifications  in  Table  II  include 
these  data,  but  the  total  measured  entry  and  the  tabulated 
momentum  spectra  do  not.  The  inclusion  increases  the  statis- 
tical accuracy  of  the  quoted  rates  and  hence  (but  only 
indirectly)  the  accuracy  of  the  spectra.  The  normalization 
factor  and  the  various  rates  have  already  been  explained  In 
the  chapter  on  determination  of  rates.  Finally,  the  tabulated 
momentum  spectra  appear.  The  various  entries  are  connected 
by  the  following  formula; 

Absolute  Intensity  - [^dlus  Intervs*^  (N°™-  p*ctor> 

(Abs.  Intensity) 

~ (No.  of  Particles)^ 


TABLE  I VERTICAL  CASE 


Direction  of  Observation: 
Magnetic  Latitude . 
Altitude 


Vertical 
48°  North 
3-3  Kilometers 


FIELD  DIRECT  DATA 


Positives  635 
Negatives  252 
Infinites,  No  Sign  3 
Total  Measured  910 
Total  Single  Tracks  1142 
Running  Time  173 o 33 
C-A  Counts  A 7co 


hrs 


Total  Frames 
Average  Yield  of 
Single  Tracks 
C-A  Rate  per  Hour 
Single  Track  Rate 
Normalization 
Factor 


2792 

41* 

27.5  + 0.4 
11.3  + 0.4 

0.0493 


Positives 


Letter  Number  of 
Group  Particles 


0- A 
A-E 
E-I 

1- L 
L-N 
N-P 
P-R 
R-T 
T-V 
V-X 
X- 


11 

57 

89 

92 

42 

69 

83 

82 

61 

47 

22 


Radius 

Average 

Interval 

Momentum 

meters 

Bev/c 

0.316 

0.048 

O.345 

0.120 

0.429 

0.213 

0.470 

0.322 

O.38O 

0.425 

0.480 

0.530 

0 0660 

0.665 

O.98O 

O.865 

1.790 

1.200 

4.150 

1.920 

Absolute  Intensity- 
particles  per  Bev/c 
cm*  sterad  hour 

1.7  + 0.6 

8.2  + l.l 
10.3  + 1.1 

9.7  + 1.0 
5.5  + 0.9 
7.1+  0.9 

6.2  + 0.7 
4.1  + 0.5 

1.7  + 0.25 
0.6  + 0.09 


0- A 
A-E 
E-I 

1- L 
L-N 
N-P 
P-R 
R-T 
T-V 
V-X 
X- 


11 

47 

88 

71 

13 

2 

3 

4 
4 

4 

5 


Negatives 


O.312 

O.330 

O.390 

0.420 

O.320 

0.390 

O.510 

O.710 

1.140 

2.150 

13.330 


0.048 

0.116 

0.200 

0.300 

0.390 

0.475 

0,585 

0.730 

0.955 

1-350 


1.7  + 0.6 
7.1  + 1.0 
11.2  + 1.1 
8.4  + 1.0 
2.0  + 0.6 
0.25+  0.1 
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TABLE  I (continued) 


FIELD  REVERSE  DATA 


Positives  1,26 

Negatives  148 

Infinite,  No  Sign  ? 

Total  Measured  477 

Total  Single  Tracks  637 

Total  Frames  1404 

Running  Time  82.64  hrs 

C-A  Counts  2404 


Average  Yield  of 
Single  Tracks 
C-A  Rate  per  Hour 
Single  Track  Rate 
Normalization 
Factor 


45% 

29.1  + 0.6 

13.1  + 0.9 

0,110 


Positives 


Letter 

Number  of 

Radius 

Interval 

Group 

Particles 

meters 

0-A 

9 

0.312 

A-E 

23 

0.330 

E-I 

44 

0.390 

X -L 

55 

0.420 

L-N 

29 

0.320 

N-P 

31 

O.390 

P-R 

31 

0.510 

R-T 

43 

O.710 

T-V 

31 

1.140 

V-X 

22 

2.150 

X- 

8 

13.330 

Average  Absolute  Intensity 
Momentum  particles  per  Bev/c 
Bev/c  cm  sterad  hour 


0.048 

0.116 

0.200 

0.300 

0.390 

0.475 

0.585 

0.730 

0.955 

1-350 


3.2  + 1,0 

7.7  + 1.6 

12.4  + 1.9 

14.4  + 1.9 
10.0  + 1.9 

8.8  + 1.6 
6.7  + 1.2 
6.7  + 1.0 

3.0  + 0.5 

1.1  + 0.3 


Negatives 


0- A 
A-E 
E-I 

1- L 
L-N 
N-P 
P-R 
R-T 
T-V 
V-X 
X- 


4 

O.316 

0.048 

27 

O.345 

0.120 

0.429 

0.213 

42 

0.470 

0.322 

12 

O.38O 

0.425 

5 

0.480 

0.530 

2 

0.660 

O.665 

2 

1. 

O.98O 

O.865 

4 

1.790 

1.200 

1 

6 

4.150 

1.920 

1.4  + 0.7 
8.6  + 1.7 
11.0  + 1.7 
9.8  + 1.3 
3-5  + 1.0 
1.2  + 0.6 
0.3  + 0.2 
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TABLE  II  EAST  CASE 


Direction  of  Observation: 
Magnetic  Latitude: 
Altitude : 


45  East  of  Vertical 

48°  North 

3.3  Kilometers 


FIELD  DIRECT  DATA 


Positives 
Negatives 
Infinite,  No  Sign 
Total  Measured 
Total  Single  Tracks 
Total  Frames 
Running  Time 
C-A  Counts 


4l6 

240 

3 

533 

767 

2099 

287.85  hrs 
2894 


Average  Yield  of 
Single  Tracks 
C-A  Rate  per  Hour 
Single  Track  Rate 
Normalization 
Factor 


36.5* 

10.1  + 0.2 
3-7  + 0.2 

0.0278 


Positives 


Letter  Number  of 
Group  Particles 


Radius 

Interval 

meters 


Average 

Momentum  Absolute  Intensity 
Bev/c  particles  per  Bev/c 
cnr  sterad  hour 


0-A 

4 

0.316 

A-E 

25 

0.345 

E-I 

39 

0.429 

I-L 

51 

0.470 

L-N 

18 

0.380 

N-P 

37 

0.480 

P-R 

39 

0.660 

R-T 

46 

O.98O 

T-V 

34 

1.790 

V-X 

26 

4.150 

X- 

14 

— 

0.048 

0.35 

+ 0.2 

0.120 

2.1 

+ 0.4 

0.213 

2.5 

+ 0.4 

0.322 

3-0 

+ 0.4 

0.425 

1.3 

+ 0.3 

0.530 

2.1 

+ 0.35 

O.665 

1.6 

+ 0.3 

0.865 

1.3 

+ 0.2 

1.200 

0.5 

+ 0.09 

1.920 

0.2 

+ 0.03 

Negatives 


0-A 

5 

0.312 

A-E 

35 

0.330 

E-I 

80 

0.390 

I-L 

59 

0.420 

L-N 

6 

0.320 

N-P 

2 

0.390 

P-R 

2 

0.510 

R-T 

- 

0.710 

T-V 

- 

1.140 

V-X 

- 

2.150 

X- 

- 

13.330 

0.048 

0.45 

+ 0.2 

0.116 

2.9 

+ 0.5 

0.200 

5.7 

+ 0.7 

0.300 

3.9 

+ 0.5 

0.390 

0.5 

+ 0.2 

0.475 

0.14 

+ 0.1 

0.585 

0.1 

+ 0.07 

0.730 

- 

0.955 

- 

1.350 

- 

SWOT* 
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TABLE  II  (continued) 


FIELD  REVERSE  DATA 


Positives 
Negatives 
Infinite,  No  Sign 
Total  Measured 


iazaL  Single  Track 
Total  Frames 
Running  Time 
C-A  Counts 


404 

238 

645 

795 

2378 

379.87  hrs 
3602 


Average  Yield 
of  Single  Tracks 
C-A  Rate  per  Hour 
Single  Track  Rate 
Normalization 
Factor 


33.4# 

9.48  + 0.15 
3.06  + 0.15 

O.OI90 


Positives 


Letter 

Group 

0- A 
A-E 
E-I 

1- L 
L-N 
N-P 
P-R 
R-T 
T-V 
V-X 
X- 


Radlus 

Number  of  Interval 
Particles  meters 


5 

27 

70 

70 

48 

58 

56 

35 

18 

11 

3 


0.312 

0.330 

O.390 

0.420 

0.320 

O.390 

0.510 

0.710 

1.140 

2.150 

13.330 


Average 

Momentum 

Bev/c 


Absolute  Intensity 
particles  per  Bev/c 
cm  sterad  hour 


0.048 

0.116 

0.200 

0.300 

O.39O 

0.475 

0.585 

0.730 

0.955 

1.350 


3.4  + 0.4 

3.2  + 0.4 

2.8  + 0.4 
2.8  + 0.4 
2.1  + 0.3 
0.9  + 0.15 
0.3  + 0.07 
0.09  + 0.05 


0- A 
A-E 
E-I 

1- L 
L-N 
N-P 
P-R 
R-T 
T-V 
V-X 
X- 


3 
49 
71 
63 
26 

8 

7 

4 
2 
4 
2 


Negatives 


O.316 

0.345 

0.429 

0.470 

O.38O 

0.480 

0.660 

0.980 

1.790 

4.150 


0.048 

0.120 

0.213 

0.322 

0.425 

0.530 

0.665 

O.865 

1.200 

1.920 


0.18  + 0.1 
2.7  + 0.4 
3.1  + 0.4 

2.5  + O.35 
1-3  + 0.3 

0.3  + 0.1 
0.2  + 0.08 
0.08  + 0.04 
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TABLE  III  WEST  CASE 


Direction  of  Observation:  45°  West  of  Vertical 

Magnetic  Latitude;  48°  North 

Altitude:  3.3  Kilometers 


(All 

data  with 

field  direct) 

Positives 

775 

Average 

Yield  of 

Negatives 

336 

Single 

Tracks 

35% 

Infinite,  N0  Sign 

16 

C-A  Rate 

per  Hour 

10.9  + 0.2 

Total 

Measured 

1127 

Single  Track  Rate 

3-8  + 0.2 

Total 

Single  Tracks 

1338 

Normalization 

Total 

Frames 

3835 

Factor 

0.0135 

Running  Time 

508.70  hrs 

C-A  Counts 

5538 

Positives 

Radius 

Average 

Absolute  Intensity 

Letter 

Number  of 

Interval 

Momentum 

particles  per  Bev/c 

Group 

Particles 

meters 

Bev/c 

cm* 

sterad  hour 

0-A 

14 

0.316 

0.048 

0.6 

+ 0.16 

A-E 

54 

0.345 

0.120 

2.1 

+ 0.3 

E-I 

101 

0.429 

0.213 

3.2 

+ 0.3 

I-L 

109 

0.470 

0.322 

3.1 

+ 0.3 

L-N 

63 

O.38O 

0.425 

2.2 

+ 0.3 

N-P 

57 

0.480 

0.530 

1.6 

+ 0.2 

P-R 

90 

0.660 

O.665 

1.8 

+ 0.2 

R-T 

96 

O.98O 

O.865 

1.3 

+ 0.15 

T-V 

78 

1.790 

1.200 

0.6 

+ 0.1 

V-X 

65 

4.150 

1.920 

0.2 

+ 0.03 

X- 

31 

“ — 

— — 

- 

- 

Negatives 

i 

0-A 

12 

0.312 

0.048 

0.5 

+ 0.15 

A-E 

66 

0.330 

0.116 

2.7 

+ O.35 

E-I 

108 

O.39O 

0.200 

3.7 

+ 0.40 

I-L 

96 

0.420 

0.300 

3.1 

+ 0.30 

L-N 

15 

0.320 

0.390 

0.6 

+ 0.20 

N-P 

8 

0.390 

0.475 

0.3 

+ 0.1 

P-R 

7 

0.510 

0.585 

0.2 

+ 0.1 

R-T 

4 

0.710 

0.730 

* 

_ 

T-V 

4 

1.140 

0.955 

_ 

— 

V-X 

11 

2.150 

1.350 

— 

— 

X- 

9 

13.330 

— _ 

• 

i ■ 

i-J 
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TABLE  IV  COMBINED  VERTICAL  DATA  1950-1951 


Direction  of  Observation; 
Magnetic  Latitude: 
Altitude : 


Vertical 
48°  North 
3 °3  Kilometers 


Positives  2666 

Negatives  1052 

Infinite,  No  Sign  31 

Total  Measured  37^9 

Total  Single  Tracks  5126 
Total  Frames  10521 

Running  Time  833 . 02  hrs 

C-A  Counts  19,878 


Average  Yield  of 
Single  Tracks  48.7# 

C-A  Rate  per  Hour  24.8  + 0.17 
Single  Track  Rate  12.1  + o!25 
Normalization  ~ 

Factor  0.0129 


Positives 


Letter  Number  of 
Group  Particles 


Radius  Average 

Interval  Momentum 
meters  Bev/c 


Absolute  Intensity- 
particles  per  Bev/c 
cm2  sterad  hour 


0- A 
A-E 
E-I 

1- L 
L-N 
N-P 
P-R 
R-T 
T-V 
V-X 
X- 


52 

199 

321 

419 

I89 

252 


309 
331 
235 
193 
1 66 


0-316  0.048 
0.345  0.120 
0.429  0.213 
0.470  O.322 
0.380  0.425 
o.48o  0.530 
0.660  O.665 
0.980  O.865 
1.790  1.200 
4.150  1.920 


2.1  + 0.3 

7.4  + 0.6 

9.7  + 0.6 
11.5  + 0.6 

6.4  + 0.5 

6.8  + 0.5 

6.0  + 0.4 

4.4  + 0.3 
1.7  + 0.1 
0.6  + 0.05 


Negatives 


0- A 

49 

A-E 

189 

E-I 

332 

I-L 

286 

L-N 

62 

N-P 

25 

P-R 

17 

R-T 

18 

T-V 

16 

V-X 

21 

X- 

37 

0.316  0.048 
0.345  0.120 
0.429  0.213 
0.470  0.322 
0.380  0.425 
0.480  O.530 
0.660  O.665 
0.980  0.865 
1.790  1 . 200 
4.150  1.920 


2.0  + 0,3 

7.1  + 0.6 
10.0  + 0.6 

7.9  + 0.5 

2.1  + 0.3 
0.7  ± 0.2 
0.3  + 0.09 
0.2  + 0.06 
0.1  + 0.03 
0.07+  0.02 


CHAPTER  V 


DISCUSSION  OF  THE  RESULTS 

For  the  vertical  case,  data  of- a similar  nature  had 
been  obtained  previously  in  studying  the  form  of  the  proton 
spectrum  (7).  Nevertheless,  observations  for  this  direction 
were  repeated,  in  part  because  of  alteration  of  the  equipment 
(Chapter  II,  Section  1)  and  m part  as  a check  on  the  previous 
work,  and  because  it  was  felt  best  that  for  comparison  between 
the  various  directions  of  observation,  the  data  should  be 
taken  In  close  succession  rather  than  a year  apart..  The 
present  and  previous  spectra  for  the  vertical  case  agree  well 
within  statistics,  a fact  that  Is  regarded  as  proof  that  the 
equipment  had  undergone  no  change  in  performance  during  the 
intervening  year  . For  improved  statistics,  the  two  sets  of 
data  have  been  combined  with  the  result  shown  in  Figure  VI. 

An  explanation  of  the  difference  in  form  of  the  spectra 
for  positive  and  negative  particles  has  been  given  previously 
in  connection  with  the  determination  of  the  proton  spectrum 
(7) . It  should  be  mentioned  here  that  the  negative  distribu- 
tion contains  mu  mesons  only,  while  the  positive  distribution 
is  for  positive  mu  mesons  plus  protons.  While  the  total 
intensity  of  positive  mesons  is  known  to  exceed  that  of  the 
negative  by  some  20#,  at  momenta  above  0.5  Bev/c  where  the 
selection  method  excludes  almost  all  mesons  the  difference 


between  the  positive  and  negative  intensities  can  be  taken  as 
a measure  of  the  proton  Intensity. 

1 . The  Meson  Spectrum 

. o 

In  Figures  II  to  V the  vertical , the  45  East,  and  the 
45°  West  spectra  are  given,  along  with  the  45°  spectrum  aver- 
aged over  East  and  West  The  latter  spectrum.  Figure  V,  is 
represented  by  plotted  points  only,  the  curves  appearing  on 
the  figure  for  comparison  purposes  being  the  combined  vertical 
data  of  Figure  VI.  The  plotted  points  are  the  measured  rates 
multiplied  by  a scale  factor  of  3.2  to  allow  comparison  of  the 
spectral  shapes.  With  this  multiplication  factor  and  within 
statistics,  the  average  East-West  spectrum  and  the  vertical 
spectrum  do  not  differ.  This  Is  true  for  both  the  proton  and 
meson  radiation.  An  explanation  of  the  zenith  angle  dependence 
of  the  meson  radiation  has  been  given  in  Charbonnler ' s study 
of  this  component.  It  was  surprising  to  find  the  same  zenith 
angle  dependence  for  protons  as  for  mesons  since  the  two  types 
of  particles  are  not  thought  to  be  closely  related  genetically 
at  this  altitude.  It  should  also  be  pointed  out  that  the  East 
and  West  spectra  agree  between  themselves  within  statistics. 
This  is  shown  in  Figure  VII.  Thus,  within  the  accuracy  of 
the  present  observations,  neither  protons  nor  mesons  at  this 
altitude  show  any  East-West  asymmetry.  The  observed  agree- 
ment in  the  zenith  angle  dependence  of  protons  and  mesons  Is 
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considered  to  be  fortuitous.  Thus  while  the  meson  intensity 
decreases  by  a factor  of  3-2  from  zero  to  45  degrees  for  the 
low  momentum  mesons  observed  here,  this  ratio  does  not  hold 
at  higher  energies.  At  2 Bev/c,  the  intensity  ratio  for 
mesons  is  found  to  be  about  2 (6) . 

It  is  not  the  primary  purpose  of  this  experiment  to 
explain  the  behavior  of  mesons.  The  zenith  angle  dependence 
of  these  particles  has  been  thoroughly  discussed  by  Charbon- 
nler  (7) . One  assumes  that  essentially  all  mesons  observed 
at  lower  altitudes  are  produced  isotropically  at  great  atmos- 
pheric heights.  Taking  Into  account  the  negligible  scattering 
of  these  particles,  their  radioactive  decay  and  ionization 
energy  loss  In  traversing  the  atmosphere,  one  accounts  quan- 
titatively for  the  observed  Intensity  and  directional  distri- 
bution at  lower  altitudes.  With  these  assumptions,  a 
transformation  can  be  found  relating  the  vertical  and  45° 
spectra  at  lower  altitudes.  In  this  way,  from  the  observed 
vertical  intensity  at  3*3  km.  for  a momentum  of  0.8  Bev/c 
one  predicts  a 45°  intensity  of  6.0  + .3  particles/Bev/c  cm2 
sterad  hour,  at  0.25  Bev/c.  The  45°  Intensity  at  0.25  Bev/c 
observed  in  the  ©resent  work  is  7.19+0.33.  Agreement 
between  the  predicted  and  observed  values  is  not  considered 
too  bad,  although  outside  the  limits  of  statistical  accuracy. 
Charbonnier  did  not  observe  meson  intensities  at  momenta 
below  0.5  Bev/c.  The  predictions  of  a 45°  intensity  at  0.25 
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Bev/c  was  nevertheless  made  since  this  depends  on  his  observed 
Intensity  at  0.8  3ev/c . At  this  low  residual  momentum,  because 
of  the  increased  importance  of  scattering,  the  computation  is 
not  expected  to  have  comparable  validity  with  Charbonnier 's 
results  at  higher  momenta.  As  the  radiation  intensity 
decreases  with  increasing  zenith  angle,  the  net  effect  of 
scattering  will  be  to  increase  the  intensity  at  large  zenith 
angles  at  the  expense  of  the  Intensity  at  smaller  angles. 

This  is  in  the  direction  of  the  observed  and  predicted 
intensities . 

2.  The  Proton  Spectrum 

As  compared  to  the  meson  radiation,  the  atmospheric 
development  of  the  proton  radiation  (and  the  nucleonic  radi- 
ation in  general)  is  greatly  complicated  by  nuclear  interac- 
tion which  here  becomes  the  most  important  single  process. 
Development  of  the  nucleonic  radiation  involves  a cascade 
process  qualitatively  similar  to  that  of  the  soft  component. 

The  cascade  is  propagated  through  the  atmosphere  principally 
by  protons  and  neutrons  (22-24)  produced  in  successive  pene- 
trating showers.  Pi  mesons  originating  in  these  showers 
decay  too  quickly  to  have  high  probability  of  undergoing 
further  nuclear  interactions  while  the  mu  mesons  secondary 
to  them  do  not  undergo  appreciable  nuclear  Interaction  at 
all,  and  so  take  no  part  in  further  propagation  of  the 
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cascade.  Compared  to  the  electromagnetic  cascade,  the  inter- 

2 

action  path  lengths  here  are  in  the  order  of  100  g/cm  , 
instead  of  the  few  g/cm  characterizing  the  radiation  lengths 
of  the  electron  cascade.  The  nucleonic  cascade  had  been 
discussed  in  detail  in  a number  of  papers  by  Heltler,  Janossy 
and  Messel  (8-20).  The  treatment  there  is  phenomenological 
in  regard  to  the  form  assumed  for  the  nuclear  cross-section. 
With  this  exception  the  mathematical  treatment  is  complete 
with  respect  to  fluctuations  resulting  from  geometrical  and 
statistical  factors.  The  cross-section  is  assumed  to  be  a 
function  of  the  fractional  energy  loss  of  the  primary,  inde- 
pendent of  the  initial  energy.  The  cross-sections  employed 
do  not  include  a description  of  the  angular  distribution  of 
the  secondary  nucleons.  Thus  the  zenith  angle  effects  have 
not  been  computed,  so  that  this  treatment  cannot  be  compared 
with  the  present  experiment . 

The  constants  in  the  cross-sections  used  by  the  above 
authors  were  so  chosen  as  to  conform  to  experimental  values 
for  the  removal  path  length  of  the  nucleonic  radiation. 
Absorption  of  the  radiation  producing  penetrating  showers 
has  been  studied  by  Van  Allen,  Bernardini,  Rossi  and  others 
(25-29,33)  who  find  an  exponential  absorption  with  a coeffi- 
cient of  about  120  g/cm^  of  air.  Studies  of  the  absorption 
of  the  radiation  producing  photographic  stars  as  well  as 
direct  measurements  by  this  laboratory  of  the  relative  proton 
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intensity  between  3.3  km  and  3ea  level  gives  similar  absorp- 
tion path  lengths.  It  should  be  pointed  out  that  an  expon- 
ential absorption  of  this  path  length  cannot  hold  to  the  top 
of  the  atmosphere,  as  in  this  case  the  observed  proton 
intensity  of  3-3  km.  would  be  much  too  large  as  compared  to 
the  primary  proton  intensity  at  the  top  of  the  atmosphere. 
Instead  (and  there  Is  some  experimental  evidence  to  the 
effect)  it  is  probable  that  for  the  very  high  energy  protons 
at  the  top  of  the  atmosphere  the  effective  removal  path 
length  Is  much  larger  due  to  the  larger  number  of  "genera- 
tions" of  secondaries  of  sufficient  energy  to  continue  the 
cascade.  In  the  following  treatment  of  the  angular  spread 
of  the  proton  radiation  a constant  removal  path  length  of 
120  g/cm  is  assumed.  In  making  comparisons  with  observa- 
tions made  at  3-3  km.,  it  is  not  necessary  that  this 
assumption  hold  to  the  top  of  the  atmosphere,  but  rather 
that  it  be  good  from  heights  of  3-3  km.  to  sea  level. 
Protons  observed  in  this  investigation  lie  principally 
between  0.5  Bev/c  and  1.0  Bev/c . On  the  basis  of  Ionization 
loss  alone,  these  protons  must  have  had  momenta  greater  by 
some  2 Bev/c  at  the  top  of  the  atmosphere.  As  they  will 
have  undergone  several  nuclear  collisions,  the  primary 
energy  must  have  been  several  times  greater  than  this.  For 
protons  of  momenta  above  7*5  Bev/c  the  primary  distribution 
is  isotropic  above  the  horizon,  as  shown  by  the  Lemaitre- 
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Vallarta  theory  (35)  of  the  geomagnetic  effect.  It  can  be 
expected  then  that  tne  protons  observed  at  3°3  km.  are 
derived  from  an  Initially  isotropic  radiation.  The  radiation 
at  3-3  km.  then  can  be  considered  as  the  superposition  of  the 
radiation  due  to  equal  contributions  from  each  direction  of 
primary  Incidence.  The  observed  angular  spread  at  3*3  km. 
as  well  as  absorption  path  lengths  as  measured  for  either  the 
total  or  vertical  radiation  depend  in  part  on  geometrical 
factors  resulting  from  this  superposition  of  contributions 
from  many  primary  directions.  It  is  the  object  of  the  fol- 
lowing treatment  to  attempt  to  separate  out  the  purely 
geometrical  factors  and  so  arrive  at  an  estimate  of  the 
angular  spread  and  absorption  pat  -.  length  as  determined  by 
the  physical  processes  and  cross-sections.  It  Is  clear,  for 
example,  that  a primary  radiation  all  of  vertical  Incidence 
would  not  ^Ive  rise  to  so  great  an  angular  spread  at  3-3  km. 
as  that  observed.  Thus  the  observed  angular  spread  repre- 
sents an  upper  limit  to  that  due  to  nuclear  scattering  in 
traversing  this  absorber  thickness.  We  attempt  to  obtain  a 
closer  approximation  to  the  angular  spread  due  solely  to  the 
nuclear  scattering  through  the  following  assumptions: 

1.  The  primary  proton  Intensity  is  the  same  at  all 
directions  of  Incidence. 

2.  The  intensity  of  secondaries  at  depth  X due  to 
primaries  of  zenith  angle  (independent  of 
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azimuth)  depends  only  on  the  diagonal  distance 

p 

d In  g/cm  as  measured  along  this  zenith  direction. 

3.  This  intensity  decreases  exponentially  with  a path 

. 2 

length  of  120  g/cm  . 

4.  The  angular  spread  of  secondaries  due  to  the  primary 
radiation  in  a given  direction  is  represented  by  a 
cos"'*)  law  centered  about  this  direction  with  n 
independent  of  depth. 

The  first  assumption  is  completely  justified  for  primary 
momenta  above  7-5  Bev/c  by  the  geomagnetic  theory.  As  a 
result  of  the  following  computation  the  assumption  that  the 
primary  radiation  is  isotropic  for  angles  greater  than  45° 
is  found  unnecessary.  For  this  smaller  angle,  geomagnetic 
theory  insures  isotropy  for  momenta  as  low  as  3-3  Bev/c. 

The  second  assumption  would  be  strictly  true  if  one 
neglected  scattering.  With  scattering,  the  geometrical  situ- 
ation is  not  strictly  the  same  for  all  zenith  angles.  A 
primary  traversing  a given  path  near  the  top  of  the  atmosphere 
will  produce  secondaries  at  lower  altitudes  with  large  lateral 
displacements  from  its  original  path.  As  compared  to  vertical 
incidence,  minimum  path  lengths  for  these  secondaries  (from 
the  point  of  observation  to  the  top  of  the  atmosphere)  will 
be  increased  for  180°  of  azimuth  and  correspondingly  decreased 
for  the  other  180°.  Except  for  large  zenith  angles,  since 
average  path  lengths  are  the  same,  the  situation  will  not 
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differ  from  that  for  vertical  incidence. 

Assumption  three  conforms  to  the  large  body  of  experi- 
mental determinations  of  this  attenuation.  While,  as  has 
been  pointed  out,  this  cannot  hold  to  the  top  of  the  atmos- 
phere, for  the  present  purpose  it  need  only  be  a close 

approximation  for  atmospheric  depths  from  about  700  to  1000 
. 2 

g/cm  . Experimental  evidence  indicates  its  validity  below 
60,000  feet  (32). 

In  regard  to  assumption  four,  a cos1^*)  scattering  law 
has  no  theoretical  foundation.  Its  use  conforms  to  usual 
practice,  n being  adjusted  to  give  the  experimentally  deter- 
mined rate  of  decrease.  The  assumption  that  n is  constant 
cannot  be  true  over  large  ranges  <f  distance.  Again  this 
only  need  be  approximately  true  over  a range  of  TOO  to  1000 
g/cm  . As  n will  be  a function  of  the  proton  energies 
involved,  similarity  of  the  proton  spectra  at  700  and  1000 
g/cm  can  be  taken  as  a measure  of  the  constancy  of  n. 
Evidence  for  this  has  been  provided  by  Wilson  (21) . 

The  physical  problem  Involved  here  is  similar  to  that 
considered  in  the  Gross  transformation  (31)#  with  the  addi- 
tion of  the  effects  due  to  scattering.  The  additional 
complexity  of  the  present  situation  necessitates  resort  to 
graphical  methods. 

Symbols  to  be  used  are  defined  as  follows  (see  Figure 
VII  for  the  coordinate  system  used) . 
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I0  {B0>^o)  Differential  proton  Intensity  at  the  top 

of  the  atmosphere  at  spherical  coordinates 

2 

6a  m partlcles/Bev/c  sterad  cm 

hour..  This  Is  assumed  to  be  constant  in 
6o  and  > 


Io  Total  differential  intensity  at  the  top 

of  the  atmosphere  « /lo(0o  cf-Oo  « 
27T  lo^fo)  particles/cm2/Bev/c  hour. 


x 


Depth  of  the  point  of  observation  measured 

2 

in  path  lengths  of  120  g/cm  . 


) Portion  of  the  differential  proton 

Intensity  at  depth  x observed  at  angles 

<9*  » % contributed  by  the  primary 

radiation  incident  at  angles  0O  , *if 0 at 

the  top  of  the  atmosphere,  in  particles/ 
2 2 

Bev/c  cm  sterad  hour. 


Ix  ($,.  <Px  ) 


Differential  proton  intensity  at  depth  x 
observed  at  angles  due  t0 

contributions  from  all  primary  radiation , 
in  particles/Bev/c  cm  sterad  hour. 


CO 


Angle  between  d-O-o  andd-Q*.  See  Figure 


VIII. 


dn0 


FIGURE  VIII 

SCHEMATIC  DIAGRAM  OF  THE 
COORDINATE  SYSTEM  USED 
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The  counting  rate  of  a Geiger  counter  telescope  at 
the  top  of  the  atmosphere  will  be  proportional  to  I (&,  <f.  )dOo 
where  da0  is  the  aperture  of  che  telescope.  The  telescope  at 
depth  x will  record  Ix(#x  > We  wish  to  account  for 

the  fact  that  this  latter  rate  falls  by  a factor  of  3.2  from 

o 

the  vertical  to  45'  zenith  angle.  For  this  purpose  we  assume 
that  ^ 

I x( 8.,  <f. ,&,><„)  = kloi&c/eo)  Ces*CO 

We  quickly  find  from  Figure  VIII  that 

£^75-  6J  - Sir)  &o  S‘r1  0^  Cos  (Wo-  *fx)  ■+•  Cos  £0,5 

so  that 

« J" kS0leo^o)e  co*&o  ^G&0oCis6lK+5w%Sm&Cos(Vo-%^f^Slx 

where  k is  a normalizing  factor  and  n0  is  the  upper  hemis- 
phere, The  determination  of  k can  be  made  immediately  by 

noticing  that  f X*  {Go^o , G)ri  <?x  ) X 

„ 

•=  klo  (Go,  ?<>)  Cos  cj  clSlK  dSl0 

JSlx 

must  reduce  to  l0{6o>  % )eH2«  for  the  case  of  no  absorption, 
l.e..  Infinite  removal  path  length.  This  means  that  x is 
zero  and  the  exponent!-1  factor  is  unity.  Thus  from 

k f \fos  G0  Cos  &x, 4 -5"?  #>  5m  &y(C°6  - ‘fx  )1  d&x  * I 

- k* 


where  J2*  is  again  the  upper  hemisphere,  k can  presumably  be 
found.  The  result,  however,  evidently  depends  on  0O  • This 
is  to  be  expected,  since  for  larger  values  of  60  , the  cosn6o 
distribution  places  particles  outside  the  limits  of  integra- 
tion. However,  it  will  be  shown  immediately  below  that  the 
exponent  n is  large  enough  so  that  the  contributions  to  the 
above  integral  for  values  of  Ok  greater  than  90°  are  negli- 
gibly small  for  0o  less  than  45°,  which  is  the  region  of 
interest.  Thus  it  is  legitimate  to  set  ■ 0,  with  the 


result 
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The  determination  of  n is  independent  of  this  result. 

Involved  in  the  choice  of  the  cosna>  dependence  for  an 
assumed  scattering  distribution  Is  the  fact  that  forward  scat- 
tering predominates  and  that  scattering  in  excess  of  90°  is 
negligibly  small.  Since  cosn«o  vanishes  at  90°,  the  behavior 
of  the  forward  scattering  is  fairly  well  approximated.  How- 
ever, cosnw  is  not  zero  beyond  90°,  and  contributions  to  the 
integral  will  occur  unless  appropriate  limits  for  the  inte- 
gration are  defined.  Hence,  the  integration  will  extend  only 
ovr  that  portion  of  the  upper  hemisphere  containing 
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directional  rays  making  an  angle  less  than  or  equal  to  90° 
with  the  direction  defined  by  , Tt  will  be  shown  later 
that  this  limitation  is  unimportant.  The  physical  situation 
demands  that  we  have  axial  symmetry  about  the  vertical  (no 
deflection  by  the  earth’s  magnetic  field  is  assumed ) , and 
hence,  the  result  must  be  independent  of  This  allows 

the  choice  - 0.  Thus  the  region  of  integration  depends 
only  on  6^  y and  Is  shown  in  Figure  .'LX,  The  shape  of  the 
region  of  integration  shows  that  if  the  coordinates  <%  , 
are  retained,  the  limits  of  integration  over  ^ will  be 
functions  of  <9  To  avoid  this  complication,  a change  of 
coordinates  is  in  order.  This  is  shown  in  Figure  X.  The 
transformed  Integral  is  then  of  ne  form.  jrx(fi.<ex)dnK 

” kJZo(0O/<{o)J  cos"( [e-  -Sir}  *Sfn  ? c/\d*l 


where  o<  = 


We  have  data  on  this  quantity  for  two  directions  of 
observation,  hence  the  integration  will  be  performed  for 
just  these  two  values,  as  it  must  be  carried  out  graphically. 

The  object  of  the  computation,  of  course,  is  to 
find  n.  The  values  of  I x(0  ) and  hav«  been 

obtained  for  n « 5,  6 and  7 at  x *=  5-4  corresponding  to  an 
atmospheric  depth  of  650  g/cm^  » 3-3  km.  altitude.  These 
are  shown  in  Table  V together  with  the  ratios  •*  I x(Of  ^ ) 


» 

J 


i*»*«i-.«**  - 


/Ix(x.  ) whl^h  are  the  measured  quantities  under  discus- 
sion., It  is  seen  that  ^ is  « fairly  sensitive  function  of 
np  and  that  by  interrogating,  n can  be  assigned  the  value  6.5. 
Also  Included  In  the  table  is  the  value  of  for  n equal  to 
M infinity,"  i.e.,  assuring  no  scattering  whatsoever  and  that 
is  simply  given  by  ^ where  x'  is  the  dis- 

tance to  the  top  of  the  atmosphere  along  the  diagonal  ray  at 
4y°  .serith  angle  to  the  vertical. 

The  value  of  n so  determined  is  dependent  on  the  choice 
of  path  length.  The  sensitivity  of  the  dependence  is  not 

great.  The  figure  for  the  path  length  as  determined  by  this  lab- 

/•* 

oratory  is  125  + 8 g /cm*.  Other  values  (29)  are  more  closely 

2 

centered  about  120  g/cm  „ Values  from  nuclear  plate  data  (26) 

2 

are  somewhat,  higher  (about  135  g/cm  ),  but  refer  to  particles 

in  much  broader  momentum  bands.  Within  the  region  L - 100  to 
2 

140  g/cm  , n varies  from  6.0  to  8.5.  Within  the  region  L » 

125  + 8 g/cm  , and  within  the  statistical  accuracy  of  the 
experimental  value  of  , n has  the  value  6,5  + 0„7. 

The  value  ■ 3,2  is  thus  shown  to  be  the  result  of 

6.5 

an  effective  scattering  distribution  of  cos  (b  . This  is  a 
fairly  sharp  distribution,  falling  off  by  a factor  of  nine 
from  zero  to  45  degrees  zenith  angle.  It  Is  of  no  small 
interest  then  that  the  use  of  such  a distribution  has  the 
effect  of  reducing  the  value  of  p by  a factor  greater  than 
three.  The  effect  of  scattering,  however,  must  be  included, 
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since  the  observed  zenith  angle  dependence  could  be  obtained 
purely  from  attenuation  considerations  only  for  path  lengths 
in  excess  of  300  g/cm2. 

The  behavior  of  the  integrand  of  the  above  integral 
shows  that  the  contribution  to  the  vertical  radiation  at 
depth  x from  the  primary  radiation  incident  at  angles  greater 
than  45  degrees  is  entirely  negligible  (less  than  0 .1%  for 
all  values  of  n greater  than  5),  and  contributions  from  the 
same  source  to  the  radiation  observed  at  45°  at  depth  x is 
quite  small  (less  than  5#  for  n greater  than  5),  primarily 
as  a result  of  the  exponential  factor.  This  provides  the 
justification  for  the  assumed  isotropy  of  the  primary  radia- 
tion, since  lack  of  isotropy  at  zenith  angles  greater  than 
45°  will  have  no  effect  on  the  computed  values  of  Ix(0  , *€ x ) 
and  Ix(^,  <fx)  . Hence,  we  need  justify  an  assumed  isotropy 
only  at  angles  less  than  45°  within  which  protons  are  iso- 
tropic at  momenta  as  low  as  3-3  Bev/c,  according  to  the 
Lemaltre-Vallarta  theory. 

Knowledge  of  the  normalization  factor  allows  one  to 
compute  an  effective  path  length  for  removal  which  takes  into 
account  the  scattering  factor.  Thus  the  computed  ratio 


Io(0  . <fo) 

I*(0.  <?») 


(both  in  the  vertical  direction)  represents  the  actual  measured 
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ratio  of  intensities  at  650  and  zero  g/cm  depths  respec- 
tively, and  when  equated  to  exp(-650/L)  yields  the  measured 
value  of  L.  These  ratios  are  exhibited  in  Table  V for  n **-  5, 
6,7  and  infinity,  as  well  as  the  associated  values  for  L. 

These  effective  values  of  L are  all  somewhat  smaller  than 
2 

120  g/cm  , as  is  expected,  since  more  particles  are  scattered 
out  of  the  vertical  beam  than  into  it,  thus  increasing  the 
apparent  attenuation  of  the  vertical  beam. 


TABLE  V 

INTENSITY  RATIOS  AND  EFFECTIVE  PATH  LENGTHS 


n 

<0>V 

_ lx(0,  V*) 

"Effective" 
Path  Length 

5 

2.28 

x 10~3 

*0 

.857  X 10-3  I0 

2.66 

IO6.5  gm/cm 

6 

2.46 

x 10~3 

I0 

.840  x 10~3  I0 

2.93 

108 

7 

2.65 

x 10~3 

I0 

.755  x 10"3Io 

3.51 

109  " 

Inf. 

4.53 

-3 

x 10  J 

I0 

.465  x 10_3Io 

9.70 

120  " 

As  mentioned  before,  the  exponential  factor  is  respon- 
sible for  the  sharp  cutoff  at  angles  in  the  vicinity  of  30-45 
degrees.  As  one  goes  to  greater  and  greater  depths,  i.e.,  to 
larger  values  of  x,  the  sharpness  of  the  cutoff  will  be 
Increased  somewhat  (in  addition  to  the  larger  attenuation) 
so  that  at  sea  level  the  result  of  a computation  of  the  ratio 
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of  the  vertical  intensity  to  the  45°  intensity  is  about  3.5. 
There  is  no  experimental  verification  of  this  surmise,.  The 
extreme  difficulty  of  picking  up  any  protons  at  all  at  sea 
level  (about  1#  of  the  total  penetrating  component)  indicates 
that  a statistically  meaningful  determination  of  this  ratio 
would  require  a long  and  careful  but  not  impractical  measure- 
ment . 

As  mentioned  in  the  introduction,  the  above  result 
places  an  upper  limit  on  the  average  deflection  suffered  by 
a proton  in  a nuclear  collision.  It  was  found  that  the  Intro- 
duction of  a co8^'(b  scattering  distribution  yielded  a satis- 
factory result  for  the  gross  zenith  angle  dependence.  A 
proton  traversing  the  distance  from  the  top  of  the  atmosphere 
to  3-3  km-  will  undergo  on  the  average  10  collisions.  This 
corresponds  to  the  geometrical  cross-section  rather  than  the 
removal  path  length  used  above.  Its  use  has  been  justified 
by  Heitler  and  Janossy  (14).  Such  protons  are  extremely  rare, 
of  course,  since  they  have  a path  length  for  removal  of  only 

twice  the  distance  corresponding  to  the  geometrical  cross- 

« 

section.  But  it  makes  no  difference  whether  one  considers  a 
single  proton  which  retains  its  identity  throughout  the  ten 
collisions  or  whether  one  asks  for  the  total  angle  that  the 
final  descendant  makes  with  the  original  primary  after  ten 
"generations,”  the  latter  situation  being  far  more  likely.  It 
is  permissible,  therefore,  to  consider  a single  proton. 
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If  we  assume  that  each  proton  suffers  exactly  ten 
collisions,  and  that  after  each  collision  the  new  direction 
differs  from  the  old  by  an  angle  of  exactly  + oC , a unique 
distribution  after  the  tenth  encounter  is  determined,  which 
depends  only  on  ot  . Since  oos^'th  was  introduced  artificially 
there  is  no  Justification  for  asking  that  this  statistical 
distribution  coincide  with  it.  One  can,  however,  compare  the 
integrals  of  the  two  distributions,  i.e.  ask  for  the  angle 
for  which  the  particle,  after  ten  encounters,  suffers  a total 
angular  deflection  of  less  than,  say,  twenty  degrees,  with  a 
certain  probability  computable  from  the  cos ^‘cb  distribution. 
This  problem  In  three  dimensions,  despite  Its  similarity  to 
the  "random  walk"  problems  of  mathematical  statistics,  and 
its  simplifying  assumptions,  is  of  enormous  complexity,  an 
unfortunate  characteristic  of  plural  scattering  problems,, 

The  two-dimensional  problem,  corresponding  to  the 
projection  of  the  three  dimensional  scattering  distribution 
on  a vertical  plane,  is  simpler,  and  actually  corresponds  to 
what  is  observed  in  this  experiment.  This  two  dimensional 
distribution  predicts  that  about  60%  of  the  particles  make 
angles  of  less  than  20  degrees  with  the  primary  direction. 

The  problem  then  is  to  adjust  0 i so  that  the  statistical 
distribution  indicates  the  same  thing.  The  choice  of  20 
degrees  as  the  matching  point  is  convenient,  for  it  is  close 
to  the  half  value  point  of  the  distribtuion.  The  percentage 
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for  a given  matching  angle  is  given  by  another  graphical 

integration,  hence  it  is  not  convenient  to  match  at  the  exact 

angle  within  which  exactly  half  the  particles  are  scattered, 

as  its  determination  is  unnecessarily  laborious.  Angles 

giving  percentages  far  removed  from  50#  are  not  expected  to 

give  the  most  accurate  matching  since  at  small  angles  the 

6.5 

projected  cos  k)  distribution  is  too  slowly  varying,  while 
at  large  angles  the  statistical  distribution  is  too  crude. 

In  a single  scattering,  it  is  assumed  that  a particle 
is  deflected  through  an  angle  of  + and  that  the  ten 

10 

encounters  are  statistically  independent.  Then  there  are  2 
possible  trajectories,  corresponding  to  six  distinct  total 
angular  deflections.  These  deflections,  along  with  their 
relative  probabilities,  are  enumerated  In  Table  VI.  One  finds 
that  a total  scattering  angle  of  about  4oc  corresponds  to  a 
probability  of  about  60#.  Hence  from  4oc  « 20°,  one  gets 
oC  - 5°. 

A check  of  this  figure  may  be  had  by  noticing  that  a 
proton's  behavior  is  similar  to  a simple  diffusion  process. 

The  root-mean-square  total  deviation  of  a proton  undergoing 
N encounters  with  a scattering  angle  of  OC  per  encounter  is 
given  by 

a = Ktf>  = T~n~ 

— o 

With  A roughly  equal  to  20  as  before,  and  with  ten  encounters. 
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oc  = 20/Xvo  — 

which  compares  with  the  previous  result,  but  the  smallness 
of  N limits  the  accuracy  of  ot  . 

Thus  a model  of  the  scattering  process  has  been  con- 
structed which  gives  a result  roughly  approximating  the 
observed  distribution.  Present  data  do  not  warrant  a more 
detailed  analysis. 


TABLE  VI 

SIMPLE  TWO-DIMENSIONAL  SCATTERING  PROBABILITIES 


e 


Relative  Probability  Probability  That  a Particle 
of  Scattering  Be  Scattered  Through  an 

at  Angle  $ Angle  Less  Than  0 


10 

8 

6 

4 

2 

0 


0 

2 

20 

90 

240 

420 

252 


1.0000 

0.9981 

0.9785 

O.89O6 

O.6563 

0.2461 

0.0000 


Total  - 1024  - 2 


CHAPTER  VI 


COMPARISON  WITH  THE  RESULTS  OF  OTHER  AUTHORS 

The  only  published  data  at  all  concerned  with  the 
problem  treated  here  is  the  work  of  Walker  (33)  and  Branch 
(36) . Walker  performed  a counter-absorber  type  of  experiment 
at  3-26  km.  elevation  and  at  sea  level  in  which  recorded 
events  are  those  due  to  charged  particles  producing  penetrat- 
ing showers  in  blocks  of  lead  absorber  of  varying  thicknesses. 
After  elaborate  precautions  to  remove  the  soft  and  the  mesonic 
components,  he  is  able  to  call  such  particles  protons,  but  the 
method  of  selection  only  vaguely  defines  accepted  momentum. 
Walker’s  requirements  for  the  lowest  energy  case  are  that  a 
charged  particle  traverse  two  trays  of  counters  separated  by 
12  inches,  enter  an  8-Inch  block  of  lead,  and  produce  a shower 
there,  whose  presence  shall  be  indicated  by  the  discharge  of 
three  or  more  counters  In  a tray  located  below  the  absorber. 
Knowledge  of  which  counters  in  the  upper  two  trays  were  dis- 
charged yields  the  zenith  angle  dependence  to  25°  In  six 
steps.  The  only  limit  that  can  be  placed  on  the  extent  of 
the  accepted  momentum  band  is  a lower  one,  which,  because  of 
the  requirement  for  shower  production,  is  certainly  very  much 
higher  than  the  1 Bev/c  maximum  cutoff  for  this  experiment . 
This  fact  is  further  substantiated  by  Walker's  ratio  of  the 
rate  at  elevation  to  that  at  sea  level.  This  ratio  is  about 


v. 
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four,  which  Is  far  removed  from  the  well-established  ratio  of 
about  twenty  as  obtained  by  this  laboratory,  by  Wilson  (21) 
and  by  Whittemore  and  Shutt  (3*0*  for  low  (less  than  1 Bev/c) 
protons . 

Walker  finds  that  his  data  at  elevation  can  best  be 
correlated  with  a cosn0  distribution  with  n equal  to  5-5  + 
0.7,  while  at  sea  level  the  exponent  is  3*^  + 1.0.  The 
increase  in  exponent  with  altitude  is  not  at  all  reasonable 
in  the  light  of  the  model  presented  here.  The  sea  level 
exponent  agrees  with  the  figure  given  in  Chapter  V,  Section 
2,  but  of  course  the  agreement  has  little  meaning  as  the  data 
do  not  overlap.  Similar  work  has  been  done  by  Brown  and  McKay 
(37). 

Branch  has  determined  the  angular  distribution  about 

the  shower  axis  of  penetrating  particles  in  extensive  air 

showers.  Such  showers,  of  course,  include,  in  addition  to 

nucleons,  many  pi  and  mu  mesons.  A conclusion  therefrom, 

relevant  to  this  problem,  is  nevertheless  possible,  since 

Branch  finds  a root  mean  square  angular  deviation  of  emergent 

particles  of  about  3 degrees.  This  figure  agrees  fairly  well 

o 

with  the  figure  5 presented  here . These  two  numbers  are  to 
be  contrasted  with  the  value  15-20  degrees  which  the  early 
theoretical  treatments  by  Messel  et  al  (8-20)  gave.  Branch's 
result,  which  this  experiment  independently  corroborates,  has 
caused  drastic  changes  in  the  nature  of  the  assumed  cross- 
section  (38,39). 


CHAPTER  VII 


SUMMARY 

The  proton  spectrum  for  particles  in  the  momentum 
band  0. 4-1.0  Bev/c  and  the  meson  spectrum  for  particles  in 
the  momentum  band  0.09-0.34  Bev/c  at  zenith  angles  45°  East* 
vertical,  and  45°  West,  have  been  obtained  at  3-3  tan*  eleva- 
tion and  magnetic  latitude  48°  North. 

The  behavior  of  the  zenith  angle  dependence  of  mu- 
mesons  has  been  found  to  be  consistent  with  a transformation 
described  by  Charbonnier  for  higher  energies,  involving  only 
the  momentum  loss  by  ionization  of  mesons  at  minimum  ioniza- 
tion and  the  well-defined  radioactive  properties  of  this 
particle . 

The  zenith  angle  dependence  of  protons,  along  with  a 
simplified  phenomenological  calculation  of  same,  is  presented 
as  new  data  on  the  nucleonic  cascade.  Theoretical  calcula- 
tions presently  found  in  the  literature  are  not  yet  well 
enough  founded  to  allow  prediction  of  these  data,  nor  have 
they  been  used  to  yield  the  differential  momentum  spectra  of 
limited  extent  as  presented  here.  It  is  found  experimentally 
that  the  gross  zenith  angle  dependence  will  fit  a cos  ' & 
distribution  at  45  degrees  East  or  West.  The  simplified 
calculation  assumes  that  a unidirectional  component  of  the 
primary  radiation  wills  (1)  suffer  exponential  absorption 


along  its  original  path  with  a path  length  for  removal  of 

p 

120  g/cm  of  air,  and  (2)  acquire  an  ultimate  distribution 
at  3-3  kw-  due  to  scattering  fitting  a cosnk>  law,  n being 
an  adjustable  parameter. 

It  is  found  that  a value  of  n of  about  6.5  yields  a 
unidirectional  distribution  which,  when  integrated  over  all 
angles,  results  in  a gross  zenith  angle  dependence  agreeing 
with  cos0'  6 . The  cos  distribution  with  n equal  to  6.5 
is  quite  sharp  and  leads  to  the  conclusion  that  the  average 
projected  scattering  angle  per  collision  is  only  about  5 
degrees . 


•.jgsaagM 
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